Quantum Dot Cavity-QED in the Presence of Strong Electron-Phonon Interactions 
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A quantum dot strongly coupled to a single high finesse optical microcavity mode constitutes a 
new fundamental system for quantum optics. Here, the effect of exciton-phonon interactions on 
reversible quantum-dot cavity coupling is analysed without making Born-Markov approximation. 
The analysis is based on techniques that have been used to study the "spin boson" Hamiltonian. 
Observability of vacuum-Rabi splitting depends on the strength and the frequency dependence of 
the spectral density function characterizing the interactions with phonons, both of which can be 
influenced by phonon confinement. 
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The fundamental system in cavity quantum electro- 
dynamics (cavity-QED) is a two-level atom interacting 
with a single cavity mode If the electric field per 
photon inside the cavity is sufBciently large, then the 
single-photon dipole coupling strength {g) between the 
atom and the cavity-mode can exceed the decoherence 
rates in the system due to cavity losses and dipole de- 
phasing: this corresponds to the strong-coupling regime 
of cavity QED whose principal signature is the vacuum- 
Rabi oscillations 0| . Cavity-QED in the strong-coupling 
regime has proved to be an invaluable tool in investi- 
gating and understanding quantum phenomena. One 
of the principal applications of cavity-QED techniques 
has been in the emerging field of quantum information 
(QI) science: a significant fraction of quantum compu- 
tation and communication schemes rely on the strong 
coupling regime of cavity-QED ||^ . Recent developments 
in semiconductor nanotechnology have shown that ex- 
citons in quantum dots (QD) constitute an alternative 
two-level system for cavity-QED applications The 
presence of electron-electron and electron-phonon inter- 
actions enriches the physics of the QD-microcavity sys- 
tem as compared to its atomic counterpart. For example, 
QD absorption/emission spectra in some cases exhibit 
sidebands and/or appreciable Stokes' shifts |^,||. These 
features signal non-perturbative electron-phonon interac- 
tions and indicate that decoherence due to phonons 
could present a fundamental limitation to QI processing 
based on quantum dot cavity-QED 0. 

In this Letter, we analyze the effects of electron-phonon 
interactions on strong electron-hole-photon coupling in 
cavity-QED, without applying Born-Markov approxima- 
tion to the electron-phonon interaction. We consider sev- 
eral forms of reservoir spectral density functions 
characterizing the interactions with different types of 
phonons. To obtain the absorption and emission spec- 
tra, we apply two related techniques that have been 
widely used to study the "spin boson" problem. The non- 
interacting blip approximation (NIB) provides an ad- 
equate solution for ohmic spectral functions {J{i^) oc uj). 



For superohmic spectral functions (J(w) oc cj",7i > 1) 
Wuerger |^ has shown that certain blip-blip interactions 
can be important, so we adopt instead a polaron opera- 
tor approach |^ that takes them into account. In both 
cases the relevant multiphonon processes are included to 
all orders and the exact solution is recovered when ei- 
ther the electron-phonon or the electron-photon coupling 
tend to zero. We find that the frequency dependence of 
J(lj) near uj — Q plays a key role in determining the 
nature of QD-cavity coupling. For superohmic reservoir 
coupling characterizing for instance exciton coupling to 
bulk acoustic phonons, we find that the strong-coupling 
regime persists even in the presence of a large Stokes' 
shift and at an appreciable temperature (T > g). How- 
ever the Rabi frequency is exponentially supressed by the 
electron-phonon coupling strength. On the other hand 
for ohmic spectral functions, vacuum-Rabi oscillations 
are only observable for Stokes' shifts smaller than a cer- 
tain critical value and for sufficiently low temperatures. 
For an ohmic continuum model this threshold value is of 
the order of the mean phonon energy {Tiojb) but if the 
phonons are strongly confined it can be much larger. For 
confined phonons in the strong-coupling regime, phonon 
satellite peaks could also exhibit splitting. 

We focus on understanding the dependence of vacuum 
Rabi oscillations on electron-phonon interactions. We 
will therefore assume a simple two-level model for the 
electronic degrees of freedom of the QD, consisting of 
the QD electronic ground state \g) and the lowest energy 
electron- hole (exciton) state |e). The starting point of 
our analysis is the Hamiltonian: 

H = huieg<Tee + hcJctt^a + hg (^(Jega + aVge) -|- 
+ (Tee J2 ^^fe {h + ^I) + E ^^f^^lK + 



(1) 



where a eg = ol and bk are annihilation operators 

for the cavity mode and the fc'^ phonon mode and fip 
a weak classical external probe that will be used for ob- 
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taining the spectra. We will also allow for Markovian 
processes, not included in to take into account cav- 
ity losses (7c) and homogeneous Lorentzian broadening 
of the zero-phonon line of the QD (7q£)). The validity 
of this model relies on the adiabatic approximation and 
the assumption that off-diagonal electron-phonon inter- 
action terms coupling |e) to exciton excited states are 
weak enough to allow the effect of the associated virtual 
transitions to be included in a dephasing contribution 
to 7q_d This assumption is well justified only for 

very small quantum dots (i? < 3 nm) and low tempera- 
tures (T < 100 K). With these approximations, we can 
describe the coupled QD-cavity system weakly excited 
by fip as a three-state system defined by projection of 
the QD-cavity Hilbert space on the subspace spanned by 
|5,ne = 0) (|0)), \g,n, = 1) (|1)) and \e,n, = 0) (|2)). Fi- 
nally, we transform the Hamiltonian to a rotating frame 
at the probe frequency. 

It is instructive to point out the relations between (|l|) 
and some well studied problems. The first three terms in 
(|l|) give the Jaynes-Cummings model in the rotating wave 
approximation |jl|]. If we project the QD-cavity system 
on any of the iV > manifolds of this model and include 
the fourth and fifth terms we obtain up to an irrelevant 
constant a spin-boson Hamiltonian |^ . Finally if we take 
(0) and drop the terms involving the cavity we get the 
independent boson model pl which has been already 
used to study phonon effects in small quantum dots . 

A convenient representation for considering super- 
ohmic electron-phonon coupling of arbitrary strength, is 
obtained by applying a canonical transformation that ex- 
actly diagonalizes the independent boson model |]^,p|,|rT| : 
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A' = e'Ae-' with s ^ (^ee^^ (bl - . (2) 



The transformed Hamiltonian reads 

H' ^Kys+Hint+HLth with: = 

k 

H'^yg = flLOaoo -I- hUcCril + fl {uJeg - A) (722 + {B)Xg 

and Hi^t^XgCg + X^Cu- (3) 



and trace over the phonon degrees of freedom 
to obtain an operator master equation for the reduced 
density matrix {p{t)) of the QD-cavity system: 



dpjt) 
dt 



E 



G„(r) 



'X,nP{t - T)e' 



(5) 



with Gg/u{t) — {£.g/u{t)£.g/u{0))- These polaron Green's 
functions, (B) and A can all be written in terms of the 
spectral function J(w) = J2k^k^i^ " ^k) The 
polaron shift A provides us with a convenient measure 
of the strength of the exciton-phonon interaction. In the 
superohmic case with n > 3, (B) = e~^l'^ where S is 
the Huang-Rhys parameter; while in the ohmic case this 
parameter is divergent and (S) = 0. 

It is lengthy but straightforward to take (||), apply per- 
turbations in and using Laplace analysis solve for the 
second order contribution to the ground state population 
(Poo (0)- '^^^ spectra can then be obtained from the 
asymptotic behavior of if) for a suitable choice of ini- 
tial condition: p(0) = (Too for absorption and /5(0) = (T22 
for emission. This choice for emission corresponds in the 
original representation to an initial condition in which 
the exciton in the QD is in thermal equilibrium with the 
phonon bath and the coupling with the cavity is suddenly 
switched on. To include the Markovian processes associ- 
ated with 7c and 7qz?, we adopt a stochastic wavefunc- 
tion (MCWF) point of view and note that the proba- 
bility amplitudes we are after coincide with the contribu- 

f2l 

tions to Pqq (t) of trajectories in which no collapse events 

occur. This implies that the spectra will still be given by 

f 2) 

the asymptotic behavior of {t) if H'^y^ is modified to 
include a non-Hermitian contribution —i^an —i-^^^<J22- 
We focus on the case where the cavity is on resonace with 
the zero phonon line (ZPL) of the QD {uj, 



eg 



and take 7c = 7qd = 7- It proves convenient to define 
Au!± — uj—ujcTg{B) and g ~ g{B), and set the frequency 
origin at uJc- For the emission spectrum we further re- 
strict ourselves to 7 = 0. In this way we obtain: 



Where we have defined the operators: 
Xg ^h[g (ct2i + (T12) + (o"20 + 0-02)] , 

Xu = ih [g {pi2 - Cr2i) + ilp (cro2 - (720 )] , 

(4) 

the mean value {B) = = (B-) and the polaron shift 
A = V, — . The Hamiltonian H',„ includes the coher- 

^-^k ujj^ sys 

ent contributions of the new interaction terms. In this po- 
laron representation we apply second order Born approx- 
imation in the residual exciton-photon-phonon coupling 



[Au;,-g^G:^{cj)Y+[^ + g^G:^'{u;)y 



(6) 



and 



G'^ (-Ai^„) G'„' (-Au;_„) 



l + e 
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+g^A'^iu;)G'^{Lo) 



(7) 



for the absorption and emission spectrum respectively, 
up to an irrelevant prefactor. Here we have defined 
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and analo- 



gous expressions for the functions G'{lu). The dissipa- 
tive parts G'^^^{iLj) and the reactive parts G'g^^{uj) are 
given respectively by the real and imaginary parts of 
Jg'^dt e^'^*Gg/„(i). The functions A'±(w) can be ex- 
pressed in terms of G'^^^{llj) and G' ^^{lo) and vanish if 
the reactive parts are set to zero Q. 

The above treatment is adequate whenever coherent 
processes, if present, are characterized by the renormal- 
ized coupling g{B). This is definitely the case for super- 
ohmic J(w) with n > 3, but for underdamped regimes 
in the ohmic case it does not hold Hence for ohmic 
J{iij) we adopt instead the influence functional formal- 
ism along the lines developed in Ref. |^ and apply the 
NIB approximation. For this case we consider only the 
absorption spectrum: the expression that we obtain co- 
incides with (||) specialized for {B) = 0. Details of this 
derivation will be discussed elsewhere []l3| . 

In physical terms the approximations we use are valid 
when the energy exhange between the QD and the cavity 
mode does not have an appreciable effect on the statis- 
tical properties of the phonons. A sufhcient condition 
for their validity is 5 ^ (5p^, where h,8ph is the smallest 
characteristic energy of J(a;) 1^,1). This condition will be 
assumed to be satisfied in all the arguments that follow. 

For superohmic J{lo) with n > 3 each of the expres- 
sions Aa;± — g^G'jJyUj) has exactly one zero {y^±) in the 
interval { — (7,5} and the functions G'^iJ^uS) and G^^^(w) 
are finite at oj = 0. In this case it is useful to write each 
of the terms 77 = — } in (^) and (0) as a "Lorentzian" 
centered at Lb±. with a frequency dependent "broadening" 
(7±(w)) and "oscillator strength" (/±(a;)). For the spec- 
tral functions and parameters we will focus on, the quan- 
titles G;/„(0)/B, G'i(c^±)/S and (a^G;/„(0))i/2 have an 
upper bound of order l/5p/i. Whenever 7±(a)±) ^ w± 
this justifies evaluating 7±(w) and /±(w) at setting 
Lb± = ±3 and dropping all contributions of the reactive 
parts. The approximation to the spectra obtained in this 
manner is analogous to the pole approximation used in 
Ref. and is valid for w ^ b^h and A{uj) (I (to)) not neg- 
ligible. We can therefore deduce that when the system 
presents weakly damped oscillations, good approxima- 
tions to the splitting, to the compound relative oscillator 
strength of the two peaks and to the broadenings are 
given respectively by 2ge~^^^, e^'^ and 7 + 2g^G'^{u)±). 
An interesting consequence of the g dependence of the 
broadenings is that the Q-value of the corresponding 
peaks is degraded if g is increased. 

The case T = 0, 7 = is of special interest. When 
T = for all J{uj), Auj- — g'^G'_{uj) has exactly one neg- 
ative zero (li-). On the other hand at zero temperature 
G'gi^{ijj) are always identically zero for lo < Q. These 
facts imply that in the limit T ^ 0, 7 ^ there is a 
lowest energy resonance at with infinite lifetime. For 
underdamped regimes at sufficiently low T and small 7 
this leads to an asymmetry between the widths of the 
two resonances (Fig. 1, 2). For the superohmic case with 



n > 3, at T = the rate 5^G"(w+) is dominated by 
one phonon processes This together with the as- 

sumptions limij^o J(w)/w" ^ ^/^ph ^^'^ '5'|t=o ~ A/cjf, 
(which are satisfied for all the J(w) we consider), allows 
us to deduce that j+/{u)+ — ui-) < {g / 5phY^~^ Se~^ im- 
plying that for T = 0, 7 = there is vacuum-Rabi split- 
ting (VRS) for aU A. 

We have computed the spectra for several specific J(w) 
when A cjf, and T < tOb- An ohmic reservoir with 
exponential cutoff (e""/"'') has been extensively studied 
in the context of the "spin-boson" problem. We find that 
at T = and for an ideal cavity, this model presents VRS 
when A/cjf, < 1 and rapid convergence to the polaron 
spectrum for A/wf, > 2 These results are consistent 
with Ref. Q , the behavior for A/w;, > 2 corresponding to 
the localization transition {a in Ref. |^ is equal to 2^)- 
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FIG. 1. Calculated spectra for a bulk phonon model with 
n = 3, A/ujb — 2, T/uib ~ 0.1 and 7 = 0. The frequency 
origin is set at the zero phonon line (ZPL). Top shows the 
polaron spectra {g = 0) and bottom the QD-cavity spectra 
for g/cdb = 5 X 10~^. The ZPL of the polaron spectra is given 
by a 5-function with relative oscillator strength = 0.325. 
The insets expand a neighborhood of = (for absorption 
the asymmetry between the two maxima is explained by the 
difference between the broadenings whilst for emission there 
is also a thermal redistribution of excitation). 

In the case of QD coupling to bulk acoustical phonons 
J{w) turns out to be superohmic. This can be deduced 
for both deformation potential and piezoelectric coupling 
|l0| and in both cases for a generic QD one finds n = 3. 
We have studied this case using a J{uj) with a power 
law cutoff given by setting J(wb)/2 = uib in equation (|^) 
presented below. It can be reasonably argued that for 
the parameter range studied g'^A'j.{iLj) <C 1 [|l3|, thus we 
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neglect these contributions. Fig. 1 shows typical spectra 
with the cavity {g ^ 0) and without it {g = 0). The 
persistence of VRS for appreciable Stoke's shift and tem- 
perature (T > g) that we find is consistent with results 
for the "spin-boson" problem |^,^. For both of these 
continuum models, ohmic and supcrohmic, Sph ~ i-^b- 

In many of the QD systems of interest strong clectron- 
phonon interactions are associated with phonon confine- 
ment ||]. To model these cases we consider coupling to 
a single phonon mode at uit which is itself broadened 
by a weak linear (position-position) interaction with a 
short memory reservoir characterized by a spectral den- 
sity J(w) with ultraviolet cutoff w*. The confined mode 
can be acoustical or optical in origin. As the interaction 
term is linear there exists an exact transformation to a 
new set of modes in which the Hamiltonian assumes once 
again the form given in (|l|) and J{lij) is given as function 
of J{i-u). The resulting J (to) inherits its ohmic or supcr- 
ohmic character from {J{uj) oc w"). For many of 
the systems of interest we can expect J{LUb) ^ and 
LUb i-^*- These conditions allow us to deduce the fol- 
lowing approximate expression for J{uj) IIT^: 



J{uj) oc 



±L + 1 



{UJ - UJb) + — ^ 



(8) 



If we assume the reservoir to be three dimensional it can 
be argued that only odd n is expected. On the other 
hand for the spectral function Sph — J{oJb)- Hence 
we have considered this model for n — 1 and n — 3 
with g <^ J{ujb)- Fig. 2 shows representative spectra for 
this model. For systems of interest one would expect: 
5 ^ 0.1 meV 0, ^ 1 mcV for acoustical phonons and 
Wb'^lOmeV for optical phonons 

An important effect of confinement on the ohmic case 
is to raise the value of critical A at which VRS is lost. 
Another striking feature is the difference between the 
ohmic and supcrohmic case as T is increased. Whereas 
for ohmic J(w) VRS is lost above a critical value of T/uJb 
that can be much smaller than unity, for supcrohmic J{lo) 
it persists for all T < uJb when A < Wb. It should be noted 
that even though is quantitatively well justified only 
for a strongly confined phonon mode, it can also be used 
as a qualitative model for an unconfined optical phonon 
band that includes optical phonon decay. In this case 
J{ujb) would be given by the bandwidth. If one considers 
the limit J{ujb) — > irrespective of n in expression (|^), 
J{lo) goes to a delta-function. For this model we find 
that the polaron operator approach constitutes a valid 
approximation and predicts for all T and A, that not 
only the ZPL, but also the sidebands present VRS ap- 
proximately given by 2ge~^^^ p^ ]. 

In summary, we have presented a formalism that al- 
lows us to analyze basic cavity-QED effects for QD- 
microcavity systems with strong exciton-phonon inter- 
actions. For supcrohmic spectral functions, the princi- 
pal role of exciton-phonon interaction is the suppression 



of the effective QD-cavity coupling strength. For ohmic 
spectral functions, non-Mar kovian decoherence due to 
phonon coupling limits the parameter range over which 
vacuum Rabi oscillations can be observed. 

We thank D. Loss for bringing to our attention Ref. 
I. W. would also like to thank B. Gayral for helpful dis- 
cussions. This work was supported by a David and Lucile 
Packard Fellowship. 
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FIG. 2. Calculated QD-cavity absorption spectra for a 
confined phonon model with A/cjb — 3, j/cob = 10"'*, 
J{uJb) — 6 X 10^^ and g/uJb = 3 x 10~^. The frequency 
origin is set at the zero phonon line and its neighborhood 
is expanded in the insets. Top shows n = 1 at T = (solid 
lines) and at T/uJb = 5 x 10~^ (dotted lines). For the latter 
the feature at w = is the remanent of an electromagnetically 
induced transparency (EIT) expected when 7 = 0. Bottom 
shows n = 3 at T/ujb = 5 x 10^^. The intensities for different 
graphs are not to scale. 
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